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STRESS  CONCENTRATIONS  IN  FILAMENT-STIFFENED 
SHEETS  OF  FINITE  LENGTH 

By  W.  B.  Fichter 
Langley  Research  Center 

SUMMARY 

A  simple  model  of  filamentary  composite  material  is  employed  to  investigate  stress 
concentrations  in  a  filament- stiffened  sheet  of  finite  length.  The  model  is  composed  of  a 
single  layer  of  parallel,  tension- carrying  filaments  embedded  in  a  shear- carrying  matrix. 
The  sheet  is  of  finite  length  in  the  filament  direction  and  of  infinite  length  normal  to  the 
filament  direction.  Filament  stress- concentration  factors  are  calculated  as  functions  of 
the  number  of  broken  filaments  and  a  length- stiffness  parameter  for  the  cases  of  uniform 
normal  edge  load  and  uniform  normal  edge  displacement.  In  the  uniform-edge-load  case, 
the  stress- concentration  factors  are  found  to  increase  with  decreasing  filament  length. 

The  opposite  effect  is  noted  in  the  uniform- edge-displacement  case  where,  in  addition,  the 
stress-concentration  factor  is  found  to  have  an  upper  limit  which  is  fixed  by  the  value  of 
the  length- stiffness  parameter. 


INTRODUCTION 

In  the  design  and  application  of  aerospace  structures  the  problem  of  stress  concen¬ 
tration  is  one  of  continuing  importance.  In  recent  years  conventional  engineering  metals 
have  been  replaced  in  some  weight- critical  applications  by  filamentary  composite  mate¬ 
rials.  These  materials  behave  somewhat  differently  from  conventional  metals  when 
weakened  by  flaws  and,  hence,  different  analytical  models  must  be  used  to  investigate 
their  behavior. 

One  simplified  model  for  filamentary-composite-material  behavior  which  has  been 
found  useful  is  the  "stringer- sheet"  model  employed  in  reference  1  for  investigating  sta¬ 
tic  and  dynamic  stress  concentrations  around  broken  filaments.  This  model  consists  of  a 
sheet  of  parallel,  equally  spaced,  tension- carrying  filaments  embedded  in  a  shear¬ 
carrying  matrix.  In  references  2  and  3  the  model  has  been  used  to  investigate  the  inter¬ 
action  of  neighboring  groups  of  broken  filaments,  the  matrix  shear  forces,  and  the 
recovery  of  load  by  broken  filaments.  In  reference  4  a  similar  model  has  been  employed 
to  study  stress  concentration  in  two-dimensional  arrays  of  equally  spaced  filaments. 


In  all  of  these  studies  attention  has  been  restricted  to  bodies  of  infinite  extent.  It 
is  well  known,  however,  that  results  for  infinite  bodies  cannot  always  be  applied  indis¬ 
criminately  to  similar  problems  involving  finite  specimens.  It  is  desirable,  therefore,  to 
assess  the  influence  of  specimen  size  on  stress  concentration.  In  the  present  paper  the 
stringer- sheet  model  is  employed  to  investigate  static  stress  concentration  due  to  broken 
filaments  in  a  filament- stiffened  sheet  of  infinite  length  normal  to  the  filament  direction 
and  of  finite  length  in  the  filament  direction.  Stress- concentration  factors  are  found  for 
a  central  cut  across  various  numbers  of  adjacent  filaments  and  for  a  wide  range  of  values 
of  a  nondimensional  length- stiffness  parameter. 

SYMBOLS 


d 

EA 

Gh 

Kr 

k 

l 

M 

Mn 

m,n,r 

N, 

P. 

P 

Pn 


n 


n 


filament  spacing 

extensional  stiffness  of  a  filament 
matrix  shear  stiffness  per  unit  length 
stress- concentration  factor  for  r  broken  filaments 
length- stiffness  parameter,  l 
half-length  of  filament  sheet 
transformed  influence  functions 

displacement  of  nth  filament  for  influence-function  solution 
indexes 

load  in  nth  filament  for  influence- function  solution 
dimensionless  load  in  nth  filament,  pn/p 
filament  load  parameter 
load  in  nth  filament 

number  of  broken  filaments 
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EA 

Ur,  dimensionless  displacement  of  nth  filament,  —  un 

11  pZ 

un  displacement  of  nth  filament 

x  coordinate  parallel  to  filament 

r /  dimensionless  coordinate  parallel  to  filaments,  - 

6  transform  variable 

Superscript: 

°°  denotes  values  for  infinite- length  sheet 

ANALYSIS 


The  finite-length  sheet  under  investigation  is  shown  in  figure  1  along  with  the  coor¬ 
dinate  system  and  some  notation.  The  filaments  are  separated  by  a  constant  distance  d 
and  are  numbered  from  -«  to  »  from  left  to  right.  The  coordinate  along  the  fila¬ 
ments  is  x  and  the  displacement  of  the  nth  filament  at  location  x  is  denoted  by  un(x). 
The  force  in  the  nth  filament,  taken  to  be  positive  in  tension,  is  pn(x)  and  is  given  in 
terms  of  un  by 


pn  =  EA 


dun 

dx 


(1) 


where  EA  is  the  extensional  stiffness  of  the  filament.  The  matrix  shear  force  per  unit 
length  between  the  nth  and  (n  +  l)th  filaments  is  given  by  —  oxn+j  -  unj,  where  Gh 
is  the  shear  stiffness  of  the  matrix.  Static  equilibrium  of  an  element  of  the  nth  filament 
requires  that 


EA  lb? + ^  (Un+1  ’ 2un  +  Un_1) =  0 


In  figure  1  filaments  0  and  1  are  shown  broken  at  x  =  0.  In  general,  for  r  broken 
filaments,  let  0  g  n  g  r  -  1  denote  the  broken  filaments.  Then,  the  appropriate  bound¬ 
ary  conditions  along  the  line  containing  the  cut  (x  =  0)  are 
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pn(0)  =  0 


un(0)  =  0 


(0  s  n  5  r  -  1) 


(n  <  0  or  n  ^  r) 


Two  sets  of  boundary  conditions  at  the  edges  of  the  filamentary  sheet  are  of  interest. 
These  are  as  follows:  (I)  uniform  normal  displacement,  and  (II)  uniform  normal  load. 
In  terms  of  the  pertinent  variables,  the  two  conditions  are 

Problem  I:  un(±Z)  =  ± 


Problem  H:  pn(±£)  =  P 

(Note  that  these  two  conditions  would  be  equivalent  if  no  filaments  were  broken.)  Prob 
lems  I  and  II  are  defined  by  equations  (2)  and  (3)  along  with  equations  (4)  and  (5), 
respectively. 

Nondimensionalization 

In  order  to  cast  the  problems  in  a  more  convenient  form,  let 


Pn  =  PPn 


Un  =  fiUn 


x  =  It 7 


Then,  equation  (2)  becomes 


_^+k2(Un+1-2Un+Un.1  =0 
dp2  v  ' 

where  k  is  a  length- stiffness  parameter  defined  by  l For  the  two  problems 
considered  here  the  complete  sets  of  boundary  conditions  become 


Pn(0)  =  0 


Un(0)  =  0 


(0  ^  n  ^  r  -  1) 


(n  <  0  or  n  ?  r) 
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and 


Problem  I:  Un(±l)  =  ±1  (9a^ 

and 

Problem  II:  Pn(±l)  =  1  (9b> 

The  dimensionless  forces  and  displacements  are  related  by 

Pn(7?)=^(l?)  (10) 

Solution  of  problems  I  and  II  is  complicated  by  the  fact  that  the  boundary  conditions 
at  x  =  0  are  mixed;  that  is,  load  is  prescribed  over  one  part  of  the  boundary,  with  dis¬ 
placement  prescribed  over  the  remainder.  To  overcome  this  difficulty  the  influence- 
function  technique  introduced  in  reference  1,  and  later  used  in  references  2  and  3,  is 
again  employed.  It  is  assumed  that  the  force  and  displacement  in  the  nth  filament  have  the 
form 

OO 

=  1  +  ^  Nn-m(,l)Um(0) 

m=  -oo 
and 


The  quantities  Mn  and  Nn  are,  respectively,  the  displacement  and  load  in  the  nth  fila¬ 
ment  due  to  a  unit  displacement  of  the  zeroeth  filament  at  77  =  0.  Application  of  the  bound¬ 
ary  conditions  at  r\  =  0  (eq.  (8))  yields 
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(11) 


r-1 


Pn(r?)  =  1  +  £  Nn_m(7])Um(0) 


m=0 

r-1 


Unfa)  =  i?  +  ^  Mn-mfa)um(°) 


m=0 


since  Um(0)  =  0  for  other  values  of  m  and 

r-1 

0  =  1  +  £  Nn_m(0)Um(0)  (0  s  n  S  r  -  1)  (12) 

m=0 

which  expresses  the  condition  of  zero  load  at  the  breaks  in  the  r  filaments  identified  by 
n  =  0  to  n  =  r  -  1.  It  should  be  noted  that  because  of  symmetry  only  the  upper  half  of 
the  sheet  (0  g  r\  ^  1)  need  be  considered.  Equations  (12)  constitute  a  set  of  r  linear 
alegebraic  equations  for  the  r  unknown  Um(0).  They  can  be  solved  and  their  solution 
substituted  into  equations  (11)  to  yield  the  entire  solution.  Before  this  can  be  done,  how¬ 
ever,  the  influence  functions  for  the  particular  problem  of  interest  must  be  determined. 
Because  of  the  similarities  between  the  influence  functions  for  the  two  problems,  they  are 
determined  here  simultaneously  in  order  to  avoid  unnecessary  repetition. 


Determination  of  the  Influence  Functions 


The  influence  functions  Mn  and  Nn  are,  respectively,  the  nondimensional  dis¬ 
placement  and  load  in  the  nth  filament.  These  result  from  displacing  the  zeroeth  filament 
a  unit  amount  at  r\  =  0  while  maintaining  zero  displacement  at  77  =  0  for  all  other  fila¬ 
ments,  and  also  while  maintaining  at  the  edge  77  =  1  either  zero  displacement  (prob¬ 
lem  I)  or  zero  load  (problem  II)  in  all  filaments. 


For  17  g  0  the  problem  can  be  stated  as 


d2Mn 

d?72 


2Mn+Mn_1)  =0 


(13) 


with  the  conditions 


Mn(0)  =  1 


(n  =  0) 


Mn(0)  =  0 


(n#  0) 


/•  \ 

Mn(l) 

few 

dq 

)  =  0 

J 

(14) 


Here  and  in  the  remainder  of  this  report  the  brace  symbol  { }  means  that  the  upper  term 
corresponds  to  problem  I  (uniform  normal  edge  displacement)  and  the  lower  term  corre¬ 
sponds  to  problem  II  (uniform  normal  edge  load). 

Let 


oo 

m (n,e)  =  ^  Mn(n)e~m6 
n=  -oo 

or,  inversely, 

Mn(^)  =  ^  J  Ml (7i,e)ein6  de 
Application  of  equation  (15)  to  equation  (13)  yields 


82M 
dr} 2 


4k2  sin2  M  =  0 


Similar  treatment  of  the  boundary  conditions  in  equations  (14)  yields 
M(O,0)  =  1 


M(1,0) 


)  =  ° 


The  solutions  satisfying  equations  (17)  and  (18)  are  found 


in 


M (77,0)  =  cosh  \2kri  [sin  | 


coth  (2k  |  sin  | 


tanh  ^2k 


sin  ^ 
2 


)sinh  r2krj 


sin  | 


(15) 


(16) 


(17) 


(18) 
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Substitition  of  these  solutions  into  equation  (16)  gives 


MM  -  ±- 


f 77 

r  / 

\ 

1 

/ 

sin  § 

\ 

coth  2k 

sln|N 

)  ( 

.  6 
Sin  2 

cosh  2k vj 

- 

\ 

/ 

2 

sin  0-  ^ 
2 

\sinh  (2krj 

\ 

tanh  (2k 

\ 

/ 

L  v 

J 

J 

-7T 


ein0  d0 


and  differentiation  yields 
77 


Nn(??)  = 


-77 


sinh  2k?] 


4k n  I  sin  ® 


- 


coth  |2k 
tanh  (2k 

v_ 


.  e 

sin- 

sin  - 
2 


\cosh  |2k r\ 


.  e 

sin  2 


.  e 
sin  2 


ein0  d6 


Finally,  making  an  appropriate  change  of  variables  and  recognizing  the  symmetry  proper¬ 
ties  of  the  integrand  yield 


tt/2 


cos  2nd  sin  6 


(coth  (2k  sin  6) 
| tanh  (2k  sin  6) 


d6 


(19) 


Stress-Concentration  Factors 

The  maximum  force  occurs  at  r,  =  0  in  the  first  intact  filament  adjacent  to  the 
broken  ones.  Therefore,  for  a  cut  across  r  filaments  the  stress- concentration  factor 
Kr  is  given  by  Pr(0)  which  is  equal  to  P_i(0)  because  of  symmetry.  Thus, 

r-1 

Kr  =  Pr(0)  =  1  +  ^  Nr_m(0)Um(0)  (2°^ 

m=0 

where  values  of  Nn(0)  for  problems  I  and  II  are  obtained  from  equations  (19).  The 
values  of  Um(0)  are  obtained  from  equations  (12)  after  substitution  of  the  appropriate 

values  for  Nn(0). 

Closed-form  evaluation  of  the  influence  functions  Nn(0)  given  by  equations  (19) 
does  not  appear  feasible,  although  asymptotic  expressions  valid  for  small  values  of  n 
and  either  very  small  or  very  large  values  of  k  are  not  too  difficult  to  obtain.  There¬ 
fore,  the  influence  functions  have  been  evaluated  by  numerical  integration  on  a  digital 
computer,  with  the  asymptotic  expressions  being  employed  in  a  few  cases  merely  as  a 
check  on  the  numerical  results.  In  addition,  the  solution  of  equations  (12)  and  calcula¬ 
tion  of  the  stress- concentration  factors  from  equation  (20)  have  been  included  m  the  auto¬ 
matic  routine. 
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RESULTS  AND  DISCUSSION 


For  both  of  the  problems  considered,  stress-concentration  factors  have  been  calcu¬ 
lated  for  various  numbers  of  broken  filaments  ranging  from  1  to  100  and  for  a  wide  range 
of  values  of  k,  the  length- stiffness  parameter.  The  results  for  problem  I  are  presented 
and  discussed  first. 


Problem  I  -  Uniform  Normal  Edge  Displacement 

The  stress-concentration  factor  Kr  due  to  various  numbers  of  broken  filaments  (r) 
is  plotted  in  figure  2  as  a  function  of  k,  the  length- stiffness  parameter.  The  results 
for  small  numbers  of  broken  filaments  are  shown  in  figure  2(a),  and  for  larger  numbers  in 
figure  2(b).  Also  shown  are  the  stress- concentration  factors  for  the  infinite -length  sheet 
(see  ref.  1)  to  which  the  present  finite-length  results  are  asymptotic  from  below  for 
large  values  of  k. 

For  very  small  values  of  k,  it  is  seen  in  figure  2(a)  that  the  stress-concentration 
factor  is  only  slightly  greater  than  one  and  is  independent  of  the  number  of  broken  fila¬ 
ments.  As  k  increases,  however,  the  number  of  broken  filaments  becomes  increasingly 
important.  For  example,  for  k  =  0. 5,  no  effect  of  the  number  of  broken  filaments  is 
evident;  whereas  for  k  =  10,  different  numbers  of  broken  filaments  yield  significantly 
different  stress-concentration  factors,  as  long  as  the  number  of  broken  filaments  is  less 
than  about  20.  For  even  larger  values  of  k,  the  number  of  broken  filaments  is  of  signif¬ 
icance  over  a  still  wider  range.  This  range  is  still  limited,  however,  for  any  finite 
value  of  k. 

In  an  experimental  program  of  tensile  tests  of  filamentary  composites  with  cuts,  it 
is  likely  to  be  k,  the  length- stiffness  parameter,  rather  than  r,  the  number  of  cut  fila¬ 
ments,  which  would  be  held  constant  throughout  a  number  of  tests.  Therefore,  it  would  be 
informative  to  present  the  results  in  a  different  way.  This  is  done  in  figure  3,  where  the 
stress-concentration  factor  is  plotted  as  a  function  of  the  number  of  broken  filaments, 
for  various  values  of  k.  Although  the  results  are  valid  only  for  discrete  values  of  r, 
they  are  plotted  as  continuous  curves  for  ease  of  illustration.  It  can  be  seen  that  for  a 
given  value  of  k  the  stress-concentration  factor  ultimately  attains  a  value  which  cannot 
be  exceeded,  no  matter  how  many  filaments  are  broken.  For  example,  for  k  =  4,  the 
stress-concentration  factor  has  an  upper  limit  of  about  2.10;  and  for  k  =  10,  an  upper 
limit  of  about  3.25.  Wherever  feasible  it  would  appear  to  be  desirable  in  many  applica¬ 
tions  to  design  a  composite  structural  component  in  such  a  way  that  an  acceptable  upper 
limit  could  be  imposed  on  the  filament  stress- concentration  factor.  This  is  analogous 
to  the  "fail-safe"  philosophy  of  aircraft  design.  Within  the  limitations  of  the  model,  this 
upper  limit  could  be  imposed  by  prescribing  a  value  of  k  and  then  suitably  adjusting 
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the  relative  values  of  length,  matrix  shear  stiffness,  filament  extensional  stiffness,  and 
filament  spacing. 


Problem  II  -  Uniform  Normal  Edge  Load 

Figure  4(a)  contains  results  for  small  numbers  of  broken  filaments;  the  results 
for  larger  numbers  are  presented  in  figure  4(b).  In  contrast  to  the  results  for  problem  I, 
the  stress-concentration  factors  for  the  uniform -edge -load  problem  are  greater  than 
those  for  the  infinitely  long  sheet,  and  they  approach  them  from  above  as  k  becomes 
large.  For  very  small  values  of  k,  it  can  be  shown  that  the  stress- concentration  factor 
Kr  approaches  the  value  1  +  This  rather  severe  stress  concentration  is  largely  of 

Ci 

academic  interest,  however,  since  uniformly  loaded  edges  are  rarely  encountered  in 
practice.  In  the  event  of  their  occurrence,  however,  the  present  results  suggest  that 
modifications  should  be  made,  perhaps  by  reinforcing  the  edges  of  the  sheet,  in  order 
that  the  boundary  conditions  would  correspond  more  closely  to  those  of  problem  I. 

CONCLUDING  REMARKS 

Filament  stress-concentration  factors  have  been  calculated  for  two  problems  con¬ 
cerned  with  the  stretching  of  a  finite-length  sheet  of  filamentary  material  weakened  by 
the  presence  of  a  group  of  broken  filaments.  For  the  problem  of  uniformly  loaded  edges, 
the  effect  of  finite  length  is  manifested  by  an  increase  in  the  stress -concentration  factor 
over  that  for  the  infinite -length  sheet.  However,  in  the  problem  of  uniform  edge  dis¬ 
placement,  the  stress-concentration  factor  is  lower  for  the  finite-length  sheet  than  for 
the  infinite -length  sheet.  In  addition,  it  is  seen  that  to  each  given  value  of  the  length- 
stiffness  parameter  there  corresponds  an  upper  limit  which  the  stress- concentration 
factor  cannot  exceed,  regardless  of  the  number  of  broken  filaments.  This  limit  should 
be  of  interest  to  designers  working  within  the  constraints  of  a  maximum  allowable 
operating  stress. 

Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Hampton,  Va.,  July  17,  1970. 
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;ure  2.-  Stress-concentration  factor  for 
r  broken  filaments  as  a  function  of  thi 
length- stiffness  parameter  k.  Uniform 
normal  edge  displacement. 


Number  of  broken  filaments,  r 


Figure  3.-  Stress- concentration  factor  for  fixed  values  of 
k  as  a  function  of  the  number  of  broken  filaments. 
Uniform  normal  edge  displacement. 


(b)  10  to  100  broken  filame 

Figure  4.-  Stress- concentration  factor  for  r  bro 
of  the  length- stiffness  parameter  k.  Unif 

NASA-Langley,  1970  -  32  L-7321 


